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(57) Abstract 

A method and apparatus for calculating the mass of a composition in a fluid volume or efficiency of exchange of said composition 
with an exchange fluid, especially urea in the body of a dialysis patient. Calculations are based on total mass of urea in the body. The 
concentration ca of urea in the effluent dialysate is measured, and the total removed mass U of urea is calculated by integrating the product 
of the urea concentration ca and the dialysate flow Qa. The momentary relative efficiency of the removal (K/V) is determined essentially by 
calculating the slope of the logarithm of the concentration curve and the. momentary mass is determined therefrom. Then, the pre-treatment 
mass of urea in the body can be determined very accurately. Moreover, the momentary relative efficiency in any point is determined by 
using the removed urea U. Hie dialysis dose is calculated by integrating the momentary efficiencies. 
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TITLE 



METHOD AND DEVICE FOR CALCULATING DIALYSIS EFFICIENCY 
FIELD OF THE INVENTION 

The present invention relates to a method and device for 
calculating dialysis efficiency using values obtained from a 
urea sensor for the calculations. The calculations can also 
predict certain conditions needing intervention. 

BACKGROUND ART 

In hemodialysis it is today common to dialyse a patient 
three times per week during a time period of three to four hours 
per treatment. The object of the treatment is to give an 
adequate dose of dialysis to the patient. Such a dose of 
treatment can be defined in different ways. 

One commonly used definition uses the urea molecule as a 
marker molecule and prescribes that the dialysis clearance (K) 
divided by the distribution volume (V) of urea times total 
treatment time (t) should exceed a certain constant, for example 
Kt/V is greater than one per treatment. The weekly dialysis dose 
is then Kt/V greater than three. 

One common way of measuring Kt/V is by measuring the 
concentration of urea (c b ) in the plasma before and after the 
treatment. The ratio R = c bpost / c bpre is correlated to Kt/V. A 
number of different equations have been suggested for the 
calculation of Kt/V, such as: 

Kt/V = -ln(R - 0,03) + (4 - 3,5-R) • UF/W (1) 
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where UF = ultrafiltration volume removed in litres and W = 
postdialysis weight in kg. 

Several clinical studies have been performed evaluating 
Kt/V in which postdialysis plasma urea c bpost has been measured 
immediately after the dialysis, usually less than two minutes 
after ending the treatment. However, most patients have a 
rebound of c bpost . If an equilibrated post treatment c bpost is 
measured after for example 30 minutes, a more "true" Kt/V can be 
measured. 

The measurement of c b is not unproblematic. It is 
required that a blood sample is taken before and after the 
dialysis treatment. Such samples are then analysed by the 
hospital's laboratory. The resulting values are given with a 
substantial time delay. In this way it is not possible to adjust 
the actual treatment so that a prescribed dose is obtained. 

The post treatment sample must be taken with care, 
especially regarding timing, to avoid false values due to 
cardiopulmonary or access recirculation. Another source of error 
is rebound mentioned above. 

If an equilibrated post treatment sample should be 
taken, such sample should be taken 30 - 60 minutes after 
terminating the treatment which is not practical for the 
patient. The amount of rebound and the rate of rebound varies 
considerably from patient to patient. 

These problems have been addressed in the prior art in 
different manners. 

WO 94/08641 describes the use of a urea monitor for 
assessing the adequacy of an dialysis treatment. The urea 
monitor is connected to the dialysis effluent line and measures 
the concentration of urea in the dialysate leaving the dialyzer. 

According to this specification, it is necessary to know 
or measure the predialysis plasma urea value (c bpre ) . Such 
measurements can be made by measuring the urea concentration in 
an equilibrated sample taken before the initiation of the 
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treatment. However, such initial measurement takes time and the 
dialysis machine needs to be specially constructed to obtain 
such predialysis urea value. 

Other indicators of adequate dialysis are URR and SRI: 

URR = 1 - R = 1 - C b post / C b pre (2) 
SRI = (m urea p re - m ure a post) / m ure a pre (3) 

where m urea pre and m urea post are pre and post amounts of urea in 
the body, respectively. 

SUMMARY OF THE INVENTION 

The object of the present invention is to provide a 
method and a device for determining the efficiency of a dialysis 
treatment and monitoring delivered dose of treatment on-line. 

Another object of the invention is to provide a method 
and device for continuously monitoring the dialysis efficiency 
for adjusting the dialysis treatment on line when required, for 
example if the dialyzer is clotted. 

A further object of the invention is to provide a method 
and device for estimating the dose of dialysis delivered without 
the need for taking blood samples or requiring the dialysis 
machine to make any special adjustments such as taking an 
equilibrated predialysis plasma urea concentration. 

According to the present invention, a urea monitor is 
. used for measuring the urea concentration c d in the effluent 
dialysate from a dialyzer and determine the total removed urea 
(U) during the treatment. The measured values are used by a 
calculation computer for estimating a predialysis urea mass m 0 
and the relative efficiency K/V. Using these values, an 
indication cf the dose of dialysis can be obtained on-line, for 
example by integrating the calculated K/V over the treatment 
time. Since the predialysis and postdialysis urea masses are 
calculated, SRI can be determined. URR can be determined as well 
if the distribution volume is estimated, for example with 
Watson's formula, see equations (2) and (3) above. Also equation 
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(1) could be used since R is known. It is noted that SRI and URR 
are obtained as equilibrated values. 

According to another approach of the .invention, it is 
assumed that the relative efficiency K/V is comparatively stable 
over at least smaller time periods and decreases continuously. 
If a sudden change in efficiency is determined, it could be an 
indication of an error condition possibly requiring nurse 
intervention, such as clotting of the dialyzer, or a change of 
blood flow Q b . 

According to still another approach of the invention, 
the effective clearance of the dialyzer is determined by 
introducing a disturbance into the dialyser and analyzing the 
effluent dialysate from the dialyser in view of the disturbance. 
The disturbance can be an alteration of the conductivity of the 
15 dialysis fluid. By the analysis of the results, it is possible 
to determine the effective clearance of the dialyser. By 
combining the dialysate concentration of urea and the effective 
clearance of the dialyser, the blood concentration of urea can 
be determined without any invasional method. 3y combination with 
20 the amount of urea obtained, by the present invention, the 
distribution volume of urea can be estimated. 

The measured concentration values of urea in the 
effluent dialysate solution has a scattered appearance for many 
reasons. However, by using a special curve adaptation algorithm, 
25 it is possible to evaluate the relative efficiency K/V over 
periods where it is relatively constant in order to accurately 
determine relevant dialysis parameters. 

Further features appears from the patent claims annexed. 
The invention and further objects, advantages and 
30 features- thereof are described in more details below with 
reference to embodiments of the invention shown on the drawinas. 
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BRIEF DESCRIPTION OF THE DRAWINGS 
Fig. 1 is a schematic view of a dialysis machine intended for 
hemodialysis including a urea monitor and where the 
invention can be used. 

Fig. 2 is a schematic view similar to Fig. 1, but with the urea 
monitor integrated in the dialysis machine. 

Fig. 3 is a schematic view similar to Fig. 1 of a dialysis 
machine adapted for hemof iltration . 

Fig. 4 is a schematic view similar to Fig. 2 of a dialysis 
machine adapted for hemof iltration. 

Fig. 5 is a diagram over concentration values obtained from the 
urea monitor in the dialysis machine according to anyone 
of Figs. 1-4. 

Fig. 6 is an estimate of initial urea mass in the diagram 
according to Fig. 5. 

Fig. 7 is a diagram similar to Fig. 5 but shows a dialysis 
treatment having a problematic portion. 

Fig. 8 is an estimate of the initial urea mass in the diagram 
according to Fig. 7. 

Fig. 9 is a schematic view similar to Fig. 2 and including 
means for inducing a disturbance in the dialyser. 

Fig. 10 is a diagram similar to Fig. 5 for determining the blood 
concentration of urea. 
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DESCRIPTION OF PREFERRED EMBODIMENTS 

The present invention is intended to be used for 
estimating parameters for a dialysis treatment, such as 
hemodialysis, hemodiaf iltration or hemofiltration. It can also 
be used for some types of peritoneal dialysis. However, the 
invention is not limited to the above-mentioned treatment modes, 
but can be used also for non-medical purposes. 

Fig. 1 is a schematic diagram of a dialysis machine 
where the present invention can be practised. The dialysis 
machine provides means for replacing the renal function of a 
mammal if the renal function is impaired or completely absent. 

The blood from a patient is taken out into an 
extracorporeal circuit 2 including a filter or dialyzer 1, 
including a semipermeable membrane 3. The blood passes on one 
side of the membrane. At the other side of the membrane, a 
dialysis fluid is circulated by the dialysis machine 4. 

The dialysis fluid is usually prepared by the machine 
from one or several concentrates and water to form a dialysis 
fluid having the desired properties. Thus, the machine disclosed 
in Fig. 1 comprises a water inlet 5, two concentrate inlets 6 
and 7, and two concentrate metering pumps 8 and 9. A first main 
pump 10 propels the ready made dialysis fluid to the dialysis 
side of the dialyzer into contact with the membrane. 

A second main pump 11 passes the effluent fluid from the 
dialyzer, the inlet dialysis fluid and any ultraf iltrate removed 
from the blood via the filter, further on to an outlet 12 and to 
the drain. 

A by-pass line 13 is arranged between the first 10 and 
the second pump 11. Several valves 14,15,16 are arranged for 
controlling the flow of dialysis fluid. The valves and the pumps 
are controlled by a computer 17 as schematically shown by 
several lines in Fig. 1. Of course, the dialysis machine is 
provided with several other means as is conventional. These 
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other means are not disclosed, since they are not relevant for 
the operation of the present invention. 

The first main pump 10 is driven with a speed so that 
the dialysis fluid delivered to the dialyzer is substantially 
constant, e.g. 500 ml/min. The second main pump 11 ia driven 
with a slightly higher speed so that the effluent fluid caUed 
the dialysate, has a flow rate of e.g. 515 ml/min. This 
operation generates a pressure at the dialysate side of the 
dialyzer, which is suitable for removing 15 ml/min of 
ultrafiltrate fluid from the blood, i.e. plasma water. During a 
treatment of 4 hours, such ultrafiltration means a fluid removal 
from the patient of 3,6 litres. Of course, the dialysis machine 
is operated so that the treatment prescribed to the patien- is 
fulfilled. 

In the effluent line from the dialysis machine is placed 
a urea monitor 18, which measures the urea concentration c d in 
the effluent dialysate. The monitor can be positioned insideTthe 
dialysis machine or completely outside the dialysis machine. The 
urea monitor can be of the type disclosed in WO 96/04401. It is 
noted that this urea monitor has a conductivity sensor, so the 
conductivity of the dialysate is determined by the urea monitor 
and the urea concentration is calculated usin such conductivity 
measurements. 

The urea monitor is shown connected to the comouter 17 
of the dialysis machine. However, the monitor can * have a 
computer of its own. 

The urea sensor or the dialysis machine also includes 
means for measuring the flow rate of the effluent dialysate Q 
The computer 17 is arranged to provide concentration values c d 
as well, as values of the total mass urea U removed during the 
treatment as the integral of Q d . Cd . The concentration values 
are taken continuously so that a concentration curve c, can be 
obtained from the urea sensor as well as a mass curve U " 
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Fig. 2 discloses a similar dialysis machine as Fig. 1. 
The main difference is that the urea monitor 19 ' is placed 
between the dialyzer 1 and the second main pump 11 and before 
the outlet of the bypass line. 

Fig. 3 discloses a similar dialysis machine as Fig. 1, 
but adapted for hemof iltration or hemodiafiltration. The only 
difference is that there is included an .infusion line 20 
including an infusion pump 21. The infusion line 20 starts from 
the outlet of the first main pump 10 and ends at the blood inlet 
side of the dialyzer, for providing an infusion fluid to the 
blood before the dialyzer, called preinfusion. The urea monitor 
22 is arranged in the effluent dialysate line after the second 
pump 11. 

Finally, Fig. 4 discloses a similar dialysis machine as 
Fig. 2, but adapted for hemof iltration or hemodiafiltration and 
providing an infusion fluid to the blood after the dialyzer, 
called postinfusion. The urea monitor 23 is placed before the 
second main pump 11 and before the outlet of the bypass line. 

Fig. 5 discloses a typical urea concentration curve c d 
obtained from the urea sensor. As appears from the figure, the 
curve is very irregular and includes several dips. These dips 
are obtained when the dialysis machine is connected for 
selfcalibration, when valve 16 is opened and valves 14 and 15 
are closed. 

Fig. 6 is a plot of urea mass values calculated 
according to the method disclosed in further details below. 

Fig. 7 is a concentration curve obtained during a 
treatment having some problematic portions as also described 
closer below. 

Finally, Fig. 8 is a plot of urea mass values calculated 
according tc the method disclosed below. 

There are several approaches to urea kinetics. One 
common approach is based on the assumption that urea is 
distributed in a single body compartment, the single pool mode^ 
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It is well known that the measured urea concentration during 
treatments does not follow such a model, specially at high 
efficiency treatments. 

Another approach assumes that the urea is distributed in 
two distinct series connected body pools with a diffusive 
exchange of urea there between. This model can explain the noted 
rebound after treatment, and the more rapid removal of urea at 
the beginning of the treatment. 

Still another approach assumes that the body is divided 
in several compartments communicating with the blood with 
different time delays. 

In the present invention, a urea monitor is used to 
measure urea concentration in the effluent dialysate from a 
dialyzer. Moreover, the total amount of effluent dialysate is 
measured. Thus, it is possible to determine the urea 
concentration c d times the total dialysate flow Q d . By 
integrating the product of c d • Q d/ the total removed urea U is 
obtained. 

If it is assumed that there is no accumulation of urea 
in the body, the total amount of removed urea (U) must be equal 
to urea generated (G) in the body over a certain time period, 
for example averaged over one week. This can be used for 
calculating the nutrition status or protein catabolic rate 
(PCR) . 

According to the present invention, the urea 
concentration measured by a urea monitor in the effluent 
dialysate solution from the dialyzer, is used for determining 
parameters of the dialysis as it progresses. These parameters 
are used for assessing the dialysis treatment on-line to 
determine the efficiency, the delivered dose, pre and post total 
urea masses in the body, urea generation rate, volume of 
distribution of urea in the body (for example by taking a blood 
sample for determining the urea concentration in the blood) , and 
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still further parameters and variables as will become evident of 
the description to follow. 

These estimations according to the present invention are 
independent of any assumptions about the distribution of urea in 
the body. 

Thus, the present invention starts from the fact that no 
assumption should be made about the urea distribution in the 
body. Instead, the total amount of urea (m) in the body at each 
moment is considered. According to the invention, the definition 
of mean urea concentration is the mean concentration of urea in 
the body over a distribution volume (V) . 

Neither the distribution volume, nor the mean urea 
concentration can be measured but we can use them for 
calculations. They can, however, be measured indirectly through 
the urea concentration measured in the effluent dialysate by the 
urea monitor as explained below. 

Moreover, according to the invention, the normally used 
dialyzer clearance is replaced by a whole body clearance K 
defined as the ratio between urea mass removal rate and mean 
urea concentration c m in the body. The urea mass removal rate is 
measured by the urea monitor and is the urea concentration (c d ) 
in the effluent dialysate times the effluent dialysate flow 
(Qd) . Consequently, these definitions are: 



K = Q d • c d / c m (4) 

(5) 



m = c m • V 



Since the rate of change of the urea mass in the body is 
urea generation rate (G) minus urea mass removal rate Q d • Cd 
the following equation is valid: 



dm/dt = G - Qd • c d = G - m • K/V 



(6) 
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This is a first order differential equation with time- 
varying coefficients which can be solved by standard methods. If 
we presume that K and G are constant and that V is a linear 
function of t (V = V 0 - UF ■ t) , we find that: 



jds = \n 



m Q-T=DF 



= -ln 



m T - 



= -ln 



KT/K, 



m 0 ~ KT/V 0 -AW/Vn 



o J 



m o-T=UF 



(7) 



The assumption made for arriving at equation (7) is that 
the whole body clearance K is constant throughout the treatment. 
10 Alternatively, it can be assumed that the relative 

efficiency K/V is constant, which results in a similar equation: 



(8) 



15 



20 



25 



However, we have found that K and K/V, respectively, are 
not constant throughout the complete treatment, but are usually 
larger during an initial period of the treatment during the 
first 30 - 60 minutes and then approximately constant. 

Moreover, equations (7) and (8) require the initial 
amount of urea, m 0 , which must be measured prior to the dialysis 
treatment, for example by a blood sample or by equilibrated 
dialysate as described in WO 94/08641 and by estimating the 
distribution volume of urea in the body. 

•However, the need for taking an equilibrated blood 
sample 30 - 60 minutes after the treatment has been eliminated 
since the equation calculates the accumulated Kt/V continuously 
so that the treatment can be terminated when the desired dose 
has been achieved. 
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Another object of the present invention is to take a 
further step in using the data obtained from the urea monitor to 
obtain all data needed for estimating important parameters of 
the complete treatment. By using the idea of whole body urea 
mass and whole body clearance K, it is possible' to develop 
equations which are independent of any assumption about constant 
clearance K or constant relative efficiency K/V over the ■ 
complete treatment period. 

By combining equations (4) and (5), and by integrating 
the left portion of equation (6) we obtain: 



and 



where 



m = (Q d • Cd) / (K/V) (9) 

(10) 



m=m 0 + G • t - U 



or 



G = urea generation which is assumed constant 
U = total removed urea, which is equal to the integral 
of Qd • c d obtained from the urea monitor. 

By rearranging we obtain: 

K/V = (Q d • c d ) / (m G + G • t - U) (ii) 
m 0 + G • t = U + (Q d • c d ) / (K/V) (12) 

Equations (11) and (12) are independent of any 
assumptions of the urea distribution in the body or any 
assumption of constant K. The only assumption made is that G is 
constant. However, if G is not constant, the product G-t should 
be replaced by the integral of G over t. 

Equation (11) can be used to evaluate the time 
dependence of K/V, which can be used for some purpose to be 
explained in more details below. 
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Like in equations (7) and (8), it is necessary to obtain 
the initial amount of urea m 0 in the distribution volume V of 
the body for equations (11) and (12). 

It is possible to obtain that value by taking a blood 
sample before the initiation of the dialysis and by estimating 
the distribution volume, for example by using Watson's formula. 

Another method is to allow the dialysate to equilibrate 
with the blood at the initiation of the dialysis treatment and 
measure the dialysate concentration, whereby the dialysate 
concentration equals the blood water concentration. By 
estimating the distribution volume, for example by Watson's 
formula, the initial mass of urea m 0 can be obtained. 

By these .methods, it is possible to obtain a value of 
the initial mass m 0 of urea, which can be used in equation (11) 
for obtaining the momentary relative efficiency K/V. By 
integrating K/V over time, an estimate of the delivered dose 
Kt/V is obtained. 

On the average, Watson's formula gives a good estimate 
of the distribution volume for a population of patients. 
However, the error for a specific patient may be large and 
difficult to predict. This may result in a large error in the 
estimate of initial mass m 0 . 

On the other hand and according to the present 
invention, it has been found that the initial mass m c can be 
calculated by using only the data obtained from a urea monitor, 
which measures the urea concentration c d in the effluent 
dialysate from a dialyzer and the flow of dialysate Q d . 

We can see that the left portion of equation (12) is a 
straight line with the slope- G. The variables c d and Q d as well 
as U (dU7dt=Q d -c d ) are measured by the urea monitor. If we know 
K/V in at least two points, it is possible to calculate the 
constants m c and G. 

It is also possible to determine G from the total 
removed urea during one week, which should be equal to the 
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generation if steady state is presumed, i.e. the generation 
equals the removal. If the patient is dialysed three times per 
week, the removed urea during one such dialysis treatment can be 
used for estimating the generation, as disclosed by Garred et 
al. 

If the constant m 0 is calculated from several measured 
values, it is possible to obtain a more accurate value of m 0 by 
taking the mean value or median value or other statistical value 
from all calculated values of the constant m 0 . 

In order to continue further, it is necessary to make an 
assumption about the time dependence of K/V. 

Firstly, it is assumed that the instantaneous relative 
efficiency K/V is constant over at least some time period during 
the treatment. The dialysis curves obtained seem to validate 
that there should be at least some periods with constant K/V, at 
least in dialysis with low ultrafiltration UF. It can be assumed 
that V is a linear function of time with constant 
ultrafiltration. However, it seems that the clearance K also 
varies with time in a similar manner, at least over certain time 
periods. 

Taking the derivative of equation (12) with constant K/V 
we obtain: 

G = dU/dt + Q d • (V/K) • dc d /dt (13) 
By inserting dU/dt=Q d -c d we obtain: 

dc d K G_ 9 K 

dt V C<i ~Q d 'v (14) 



By integrating and taking the logarithm of equation (14) 
we obtain: 



ln(c~-G/Q d ) = ln(c c -G/Q d ) - Kt/V ( i 5) 
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The curve of equation (15) is a straight line with the 
slope K/V. As appears from equation (15), the dialysate 
concentration value c d has to be reduced by an offset term G/Q d 
related to the generation. 

Thus, by using equation (15) in a period where the 
slope of the curve is constant, the momentary relative 
efficiency K/V can be determined in a number of points. Thence, 
by using equation (12), m 0 can be calculated for each point and 
an averaged (mean or median) value of m 0 can be estimated. 

Secondly, it is assumed that K is constant over at least 
some time period during the treatment. The dialysis curves 
obtained seem to validate that there should be at least some 
periods with constant K. We also presume that during said time 
period the distribution volume of urea in the body V is a linear 
function of time t with a constant ultrafiltration UF: 



V = V 0 - UF • t 



(16) 
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By taking the derivative of equation (12) we obtain: 



G = dU/dt + (Q d -V/K) -dc d /dt - (Qd-c d /K) -dV/dt (17) 



25 



which results in: 



dCd/dt + c d • (K -UF)/V = (K/V) • G/Q d 



(18) 



Equation (18) can be solved: 



30 



cAO- 



C(0)- 



QA\-UFIK) 



K-UF 
UF 



(19) 
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The expression G/[Q d (1 -UF/K) ] - c k is an offset term 
for the measured dialysate concentration c d due to urea 
generation G. The dialysis concentration c d will approach this 
offset c k asymptotically for long treatment times. 

Thus: 

ln(c d -c k ) = ln(c 0 -c k )+ [ (K-UF) /V 0 ]- [ (V 0 /UF) In (V/V 0 ) ] (20) 

where 

[ (V 0 /UF)ln(V/V 0 ) ] = (V 0 /UF)ln(l-t-UF/V : ) s -t (21) 



By plotting the left member of equation (20) versus 
[ (Vo/UF) ln(V/V 0 ) ] it is possible to determine the slope 
[(K-UF)/V 0 ]. The ultrafiltration UF is constant and known. Thus, 
the relative efficiency K/V c can be determined, if UF/V 0 is 
15 estimated. By using equation (12), m 0 can be determined, if the 
urea generation (G) is known. Otherwise, the urea generation can 
be determined as the time varying component of the determined 
m 0 . 

The determination of m 0 is made only over the time 
20 period where it is found that K or K/V, respectively, are 
constant, i.e. where the measured data, fit equations (15) cr 
(20) well enough. 

As can be seen from the enclosed diagram, Fig 3, over a 
dialysate concentration curve from a typical patient, c d is very 
25 irregular due to inter alia bypass periods, cell-to-cell checks 
in the urea monitor, noise and changes in the treatment 
efficiency (K/V) for various reasons. 

It is not easy to find a curve which is the best fit for 
a certain time period where K or K/V should be constant. 
However, the Hough transform (US 3069654), used for finding' 
lines in images, is a method capable of handling such types of 
disturbances. The Hough transform looks for straight lines that 
passes through the largest number of points. So even if a large 
number of points are outside the line, this method can still 
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work. It will also find several lines,, if there are changes in 
the treatment efficiency. 

According to the preferred embodiment of the present 
invention, the following steps are performed. 
5 The urea generation rate G is estimated from patient 

data, for example removed amount of urea during a week. 

A continuous measurement is performed in the effluent 
dialysate of urea concentration c d . At the same time the urea 
concentration c d times the dialysate flow rate Q d is integrated 
10 to provide the total removed urea U. 

The measurement is started at time zero, which is 
defined as the time where the measurement of urea concentration 
exceeds a certain level during more than five minutes. 

From time zero, the urea concentration c d is plotted 
15 versus time and the total removed urea U is calculated by 
integrating the urea concentration c d times Q d . 

Then, there is a waiting period of for example 60 
minutes, where it is assumed that K or K/V may be changing. 

After the waiting period, the data of the urea 
2C concentration curve is processed by subtracting the offset term 
G/Q d from the urea concentration c d and plcrting the logarithm 
of the corrected urea concentration. Then, the curve is 
processed for. finding a portion where K/V is substantially 
constant, as discussed above in connection with equation (15). 
25 This is performed by using the Hough transform of finding a line 
which passes through the largest number of points on the 
logarithm of the corrected urea concentration curve. When a 
portion of sufficient length has been located, the slope of the 
curve is determined for calculating K/V. 
30 Then, a number of measurement values of the urea 

concentration c d is selected which are within a certain 
deviation from the line obtained by the Hough transform, for 
example within 1% from the line. For these points, the 
instantaneous mass of urea is calculated by using equation {S. 
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Finally, these instantaneous mass values are referenced to time 
zero, as defined above, by using equation (10) for obtaining 
several calculated values of m 0 . The median value of these 
calculated initial mass values are regarded as. the best estimate 
of the initial mass m 0 of the body. 

Our research has shown that the estimate of the initial 
urea mass is very accurate. Thus, the estimate of the efficiency 
and dose of the treatment based on this invention are also very 
accurate . 

When the initial urea mass m 0 has been obtained, it can 
be used in many different ways to estimate the efficiency of the 
treatment. 

The dose of the treatment can be calculated by using 
equation (11) to obtain K/V at each moment of time. Then, K/V is 
integrated over the time to obtain Kt/V. When the desired dose 
has been obtained, the dialysis treatment is terminated. 

If the distribution volume V 0 of urea in the body is 
estimated, the initial urea blood concentration c b pre in the 
body can be obtained without the need for taking any blood 
sample by dividing the estimated mass m 0 with the distribution 
volume V 0 . Since the total removed urea U is calculated 
continuously, the urea mass after the dialysis treatment is 
known and the post urea blood concentration c b post can be 
estimated since ultrafiltration is also known. The post urea 
concentration is the equilibrated urea concentration, since the 
method according to the present invention calculates on total 
urea mass in the body. 

Thus, URR or SRI can be calculated according to 
equations (2) or (3). Moreover, equation (1) can be used for 
estimating the dialysis dose Kt/V. 

The obtained determination of the initial mass m 0 can be 
used for calculating the urea distribution volume V of the 
patient, for example after completion of the dialysis session. 
First, the initial urea concentration in blocd c b pre is measured 
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before the dialysis to obtain a value of the mean urea 
concentration before dialysis, for example by a blood sample or 
equilibrated dialysate measurement. This urea concentration is 
equal to the "mean" urea concentration according to the present 
invention, since urea is initially equally distributed in the 
body of a patient. Then, the distribution volume V pre (V 0 )can be 
calculated by dividing the initial urea mass m 0 as calculated 
according to the present invention with measured blood 
concentration c b pre . Finally, the ultrafiltration volume removed 
during the treatment is subtracted to obtain V post , i.e. the 
distribution volume of urea after the dialysis treatment. This 
post distribution volume Vpo„ should be fairly constant for a 
normal patient in a steady state and can be used as an 
additional clinical parameter. 

There are many alternative methods of using the 
principles of this invention. For example, it is possible to use 
the alternative method where K is assumed constant as discussed 
in relation to equations (16) to (19) . 

Instead of using the urea concentration values which are 
within 1% from the Hough transform line, the line itself can be 
used for the calculations. Moreover, the total removed urea U 
curve can be approximated with one or several exponential curves 
using the Hough transform. 

Other types of curve adaptation algorithms can be used. 
Thus, it is possible to use the least squares method. In that 
case, it is necessary to remove the data portions where the 
dialysis machine has made some self-calibration etc. This can be 
done in an iterative manner, where a first approximation is made 
and all data portions outside a certain limit, such as 10%, are 
removed and the process is repeated. 

In certain types of dialysis machines, the dialysis 
treatment is interrupted regularly for self-calibration of the 
machine. In that mode, valve 16 is opened while valves 14 and 15 
are closed, see Fig. 1. Thus, the dialysis cf the blood ceases 
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after a short while when the dialysis solution inside the 
dialyzer has obtained equilibrium with the blood. In the urea 
concentration curves such self-calibration periods are visible 
by regular dips in the curve, see Fig. 5, where such self- 
calibration is carried out with 30 minutes intervals. After each 
such self-calibration, the dialysis starts again at a slightly 
higher level. 

In order to account for such intermittent stops in the 
dialysis, the time scale should be adjusted to remove a portion 
of the stop period, since there is obviously no dialysis at 
least during a portion of the . stop period. We have found that 
the best approximation to reality is obtained if the stop period 
is replaced by a period of 30 seconds. This is substantially 
independent of the actual length of the stop period, which can 
be anything from 35 seconds to several minutes. 

As explained above, it is possible with the present 
invention to obtain a value of the initial mass of urea in the 
body of a patient. The invention can also be utilised in other 
areas, where it is of interest to know the mass of a substance 
or composition, such as in a beer brewery. 

Other substances than urea can be monitored, such as 
creatinine, sodium, potassium, calcium, magnesium, bicarbonate, 
glucose, 15 2 -microglobuline, etc. It is also possible to monitor 
the conductivity of the plasma water or blood or the osmolarity 
thereof. It is also possible to use the principles of the 
invention in connection with gases, such as oxygen gas, nitrogen 
gas or carbon dioxide gas. 

If the invention is to be used for compositions in the 
body having some active mechanism interfering in the body, such 
as for sodium and potassium ions, such interaction should be 
taken account for. 

For sodium and potassium and some other solutes, it is 
customary to include some concentration of these ions in the 
fresh dialysis solution and therefor it is necessary to 
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calculate on the difference between initial concentration and 
final concentration in the effluent dialysate from the dialyzer. 
One approach is to replace the concentration value c d by the 
concentration difference c dout minus c dir and the mean 
concentration c m by the difference between the mean 
concentration and the initial dialysis fluid concentration, i.e. 
c m minus c din . Consequently, equation (4) at page 9 will 
essentially be replaced by: 

K = Q d • (c dout - c doin ) / (c m - c din ) (22) 

If the treatment has followed a standard treatment with 
no apparent complications, the dialysis dose can be calculated 
by using the initial and final dialysate urea concentrations 
c d pre and c dpost and using the equation: 

URR = 1 - C d post/ C d pre (23) 

If the URR calculated using the method according to the 
present invention differs substantially from the URR obtained 
with equation (23), it is an indication of problems during the 
dialysis, such as clotting of the dialyzer, which otherwise 
could have passed undetected. 

For -determining the urea mass by equations (15) or (19) 
it is assumed that the concentration follows an exponential 
curve over at least a portion of the curve. Since the removed 
urea mass U is the integral of the concentration c d multiplied 
by the dialysate flow Q d (which is constant;, it follows that 
also U is an exponential curve over at least a portion thereof. 
Consequently, it is possible to use U instead of c d for 
calculating the momentary relative efficiency. 

Alternatively, U can be used for verifying that the 
calculations using the concentration c d are correct. Thus, it 
can be assumed that. U approaches an asymptote which is: 
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Asy = m 0 + G-t - G / (K/V) (24) 

Since all these constants are obtained via the equations 
5 given above, it is easy to calculate U minus Asy to see if this 
curve is an exponential curve with the same exponent as the 
concentration curve. If this is not the case, there is probably 
some error. 

According to the invention, the initial mass of urea is 
10 determined and several clinical parameters are calculated 
therefrom. However, the blood concentration of urea cannot be 
obtained, but needs to be measured by taking a blood sample and 
analysing it later, or by equilibrated ultrafiltration before 
the dialysis treatment is started. However, it is possible to 
15 determine the effective clearance of the dialyser with a method 
where a distrubance is introduced into the dialyser and the 
resultant effect on the effluent dialysate is analysed. 

Such a method is shown in Fig. 9 which is a schematic 
view similar to Fig. 2. In the dialysis circuit is added a pump 
24 connected to the inlet of the dialyser between the valve 14 
and the dialyser 1. To the other side cf the pump 24 is 
connected a bag 25 comprising a material to be added to the 
dialysis circuit via pump 24. 

Moreover, Fig. 9 shows a pump 26 connected to the blood 
25 circuit at the inlet of the dialyser 1 for introducing a 
material comprised in a bag 27 connected tc the other side of 
the pump 26. 

Any of these devices can be used for introducing a 
disturbance to the inlet of the dialyser. It is also possible to 
produce • a disturbance by operating the concentrate pumps 8 
and/or 9. 

The disturbance is a change of a parameter of the 
dialysis fluid or the blood. The disturbance can be a change of 
the conductivity or a change of the urea concentration. It is 
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noted that the urea monitor can measure both urea concentration 
and conductivity in the effluent dialysate. if another 
measurement instrument is used, any other substance can be used 
as a disturbance, as soon as it is compatible with the body, 
such as sodium, bicarbonate etc. 

The influence by the dialyser on the disturbance is 
measured downstream of the dialyser, for example by the urea 
sensor. A portion of the disturbing material will pass the 
membrane from the dialysate to the blood cr vice versa. The 
amount passing the membrane is dependent ■ on the dialysance of 
the membrane. 

If the disturbance is a step change in the conductivity, 
produced by pumps 8,9 the dialysance of the dialyser can be 
determined according to equation (see EP 547 025, the contence 
of which is included in the present application by reference) : 

De = Qd [1 - (Cdout2-C do utl) /(Cdin2-Cciinl)] (25) 

where 

D e = effective dialysance of the dialyser 
Q d = effluent dialysate flow 

Cdouti and Cdout2 = concentration in the effluent dialysate 

Cdini and c din2 = concentration in the introduced dialysis fluid 

Indexes 1 and 2 indicates before and after the step 
change. The introduced concentration can be measured or be 
determined by the set values of the concentration pumps. 

it is also possible to determine the effective 
dialysance of the dialyser by the method disclosed in EP 658352 
where three concentrations are measured and the dialysance is 
determined as disclosed in said patent specification, the 
contence 'thereof being included in the present application by 
reference. 

An alternative method of determining the effective 
dialysance is disclosed in EP patent application No. 9715818, 
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the contence of which is included in the present application by 
reference. The dialysance is determined by the formula: 

D e = Q d x (1 - S out /S in ) (26) 

where: 

D e = effective dialysance of the dialyzer 

Q d = dialysate flow emitted from the dialyzer 

S out = integral of Qd x (cd(t) - cdO) during the disturbance in 

the flow emitted from the dialyzer 
S in = integral of Qd x (cd(t) - cdO) during the disturbance in 

the flow entered into the dialyzer 

The disturbance can be a change of conductivity or a 
change of urea concentration or any other substance that can be 
measured and is compatible with the body. 

After obtaining the effective clearance of the dialyser 
with any of the above-mentioned method, it is observed that urea 
monitor measures the concentration of urea in the effluent fluid 
continuously. Thus, the urea concentration at the start of the 
treatment can be extrapolated from the first 5 to 20 minutes of 
treatment as shown in Fig. 10. Then, the plasma water 
concentration of urea in the body at the start of the treatment 
can be determined according to the formula: 

Cpw = Qd X Cd / K e / o 7 v 

where K ] 

c pu =plasma water concentration of urea at initiation of dialysis 
Qd = effluent dialysate flow rate 

c d = concentration of urea as extrapolated to the initiation 
K e = effective clearance of the dialyser for urea 

Since the plamsa water concentration of urea can be 
calculated as indicated above and the amount of urea at the 
start of the treatment is estimated according to the present 
invention, the distribution volume V of urea in the body can be 
calculated. This distribution volume V is an important clinical 
parameter, which now can be measured with high accuracy. 
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The invention has been described in connection with 
removing a substance from the body, such, as urea. The same 
principle is valid for the addition of a substance to the body, 
such as bicarbonate, acetate or lactate, etc. 
5 The invention has been described in connection with a 

urea monitor, which measures the urea concentration in the 
dialysate continuously. It is also possible to use a measurement 
apparatus which measures the concentration intermittently, for 
example with one or a few minutes interval. 
10 In principle, the invention can also be used for 

peritoneal dialysis, where the effluent dialysate is monitored 
for a certain substance or composition. Specially at tidal 
automatic peritoneal dialysis, where the dialysate in the 
patient is partially replaced periodically, the principles of 
15 this invention could be applied. 

The invention has been described above with reference to 
the embodiments shown in the drawings. The various components 
and characteristics can however be combined in different ways 
than have been shown in the drawings and other combinations are 
included within the scope of the invention. The invention is 
only limited by the appended claims. 



2G 



WO 98/55166 



PCT/SE98/01048 



26 



CLAIMS 

1. A method for calculating the mass of a composition in 
a fluid volume, comprising: 

performing mass exchange between said volume and an 
exchange fluid flow, 

repeatedly measuring the concentration of said 
composition in the exchange fluid flow after said mass exchange 
for obtaining a concentration curve, 

characterized by 

fitting an approximation curve to at least a portion of. 
said concentration curve, where the logarithm of said 
approximation curve is a substantially straight line, 

determining a parameter of said approximation curve for 
use in calculating the mass of said composition in said volume. 

2. A method as claimed in claim 1, characterized in that 
said composition in a fluid volume is a solute in a liquid 
volume. 

3. A method as claimed in claim 2, characterized in that 
said mass exchange is performed by passing said solute and 
liquid on one side of a semipermeable membrane in a dialyzer and 
said exchange fluid is a dialysis liquid passing on the other 
side of said membrane resulting in a dialysate liquid exiting 
the dialyzer. 

.4. A method as claimed in claim 2, characterized in 
calculating the accumulated mass' of said solute in said 
dialysate liquid by measuring the flow rate Q d of said dialysate 
liquid flow and integrating the product of said flow rate Q d and 
said solute concentration difference c d over time. 
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5. A method as claimed in claim 3 or 4, characterized in 
that said dialysis liquid entering the dialyzer comprises a non- 
zero initial concentration of said composition, and that 
measuring of the concentration is replaced by measuring the 
concentration difference across the dialyzer. 

6. A method as claimed in claim 3, 4 or 5, characterized 
in that said fitting of an approximation curve involves 
subtracting a compensation term (G/Q d ) from said solute 
concentration c d to obtain a compensated concentration, taking 
the logarithm of said compensated concentration and fitting a 
straight line to said logarithm of said compensated 
concentration, where said compensation term compensates for 
generation (G) of said solute. 

7 . A method as claimed in any one of the preceding 
15 claims, characterized in that said parameter of said 

approximation curve is a slope of said substantially straight 
line. 

8. A method as claimed in any one of the preceding 
claims, characterized in that said fitting is excluded for data 
obtained during an initiation period, of for example 60 minutes. 

9. A method as claimed in any one of claims 3 to 8, 
characterized in that said flow of dialysis liquid is 
periodically interrupted for a first time period, and the time 
scale is adjusted by replacing said first time period with a 

25 replacement time period which is shorter than said first time 
period. 
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10. A method as claimed in any one of claims 3 to 3, 
characterized in that said parameter is the slope (K/V) with 
respect to time (t) of equation: 

ln(c d -G/Q d ) = ln(c 0 -G/Q d ) - Kt/V (I) 

where 

c d = dialysate concentration at time t 
G = solute generation 
Q d = dialysate flow 

c 0 = dialysate concentration at time zero 
K/V = relative dialysis efficiency 
t = time from time zero. 

11. A method as claimed in any of the preceding claims, 
characterized in that a momentary mass (m x ) is determined 
according to the equation: 

m x = (Q d • c dl ) / (K/V) i 
where (K/V) 1 is determined for exampel according to equation 
(I), and in that the momentary relative efficiency (K/V) 2 at any 
time is determined according to equation: 
(K/V) 2 = (Q d • c d2 ) / m 2 

where 

m 2 = mi - (U 2 - Ui) + G (t 2 - ti) (iv) 

and 

c d i = concentration at time t 2 
c d2 = concentration at time t 2 
Ui = accumulated mass to time t 2 
U 2 = accumulated mass to time t 2 . 

12. A method as claimed in claim 11, characterized in 
that the momentary relative efficiency (K/V) is integrated over 
time to give an estimate of the total dialysis dose Kt/V. 

13. A method as claimed in claim 11 or 12, characterized 
in that the fitting of said approximation curve is performed by 
calculating a line which passes through the greatest number of 
points in the logarithm of the concentration curve, possibly 
compensated. 
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14. A method as claimed in claim 13, characterized in 

that all values of concentration which are within a certain 
limit from said line are used for calculating the momentary 
masses m n , which are then used for calculating the initial 
masses m 0n , and in that the calculated initial masses mo„ are 
used for estimating the actual initial mass m 0 , for instance by 
taking the median or mean value of the calculated initial masses 
n\o n - 

15. A method as claimed in claim 11 or 12, characterized 
in that the determination of a momentary mass m a is performed in 
another way, for example by analysing a blood sample or by 
equilibration of dialysate with blood and determining the actual 
concentration of the composition and by estimating or measuring 
the distribution volume of the composition. 

16. A method as claimed in anyone of claims 3 - 15, 
characterized in that the distribution volume (V) of the 
composition fluid is estimated, for example by Watson's formula, 
and that the concentration of the composition in the fluid is 
determined by dividing the calculated mass with the volume. 

17. A method as claimed in anyone of claims 3 - 15, 
characterized in that the concentration of the composition In 
the fluid is measured and that the distribution volume is 
determined by dividing the calculated mass and the 
concentration. 

18. A method as claimed in anyone of claims 3-15 
characterized in that the concentration of the composition in 
the fluid is measured by introducing a disturbance into the 
dialyser and measuring the resulting effect in the efflue- 
dialysate and calculating the effective clearance of the 
dialyser from this resulting measurements, calculating the 
plasma water concentration of said substance by the formula: 

Cpw = Qd x c d / K e 

where 

c pw =plasma water concentration of urea at initiation of dialysis 
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Qd = effluent dialysate flow rate 

c d = concentration of urea as extrapolated to the initiation 
K e = effective clearance of the dialyser for urea 
and determining the distribution volume V of said substance by 
the formula: 

V - m 0 / c pw 

where m 0 is determined according to any of the preceeding claims 

19. A method as claimed in anyone of the preceding 
claims, characterized in that, after establishment of said 
approximation curve and determining a deviation of the 
concentration curve from the approximation curve, an alarm is 
emitted at deviation over or under a predetermined threshold 
level . 

20. A method as claimed in any one of claims 3 to 9, 
characterized in that said parameter is the slope [(K-UF)/V 0 ] of 
ln(c d -c k ) as a function of (V 0 /UF) In (V/V 0 ) in the equation: 



where 



ln(c d -c k ) = ln(co-c k )+ [ (K-UF) /V 0 ]- [ (V 0 /UF) In (V/V 0 ) ] 



(V) 



c d = dialysate concentration at time t 

c 0 = dialysate concentration at time zero 

Ck = G/[Q d (1 -UF/K)] 

G = solute generation 

Q d = dialysate flow 

K = whole body clearance 

V 0 = distribution volume before treatment (pre) 
UF = ultrafiltration per time 

21. A method as claimed in claim 2C, characterized in 
that a momentary mass ( mi ) is determined according to the 
equation: 

n>i = (Qd • Cdi) / (K/V)i (VI) 



where 

(K/VJ: = (K/V 0 ) / (1 - t a • UF/ V 0 ) 



(VII) 
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with (K/V 0 ) determined according to equation (V) and UF/V 0 
estimated, and in that the momentary relative efficiency (K/V) 2 
at any time is determined according to equation: 

(VIII) 



where 
and 



(K/V} 2 = (Q d • Cd2 ) / m 2 
m 2 = mi - (U 2 - Ui) + G (t 2 - ti) 



(IX) 



c d i = concentration at time ti 
c d 2 = concentration at time t 2 
Ui = accumulated mass to time ti 
U 2 = accumulated mass to time t 2 . 

22. An apparatus for calculating the mass of a 
composition in a fluid volume, comprising: 

an exchanger for performing mass exchange between said 
volume and an exchange fluid flow, 

a measuring apparatus for repeatedly measuring the 
concentration of said composition in the exchange fluid flow 
after said mass exchange for obtaining a concentration curve, 

characterized by 

firsz calculation means for fitting an approximation 
curve to at least a portion of said concentration curve, where 
the logarithm of said approximation curve is a substantially 
straight line, 

second calculation means for determining a parameter of 
said approximation curve for use in calculating the mass of said 
composition in said volume. 

23. An apparatus as claimed in claim 22, characterized 
in that said composition in a fluid volume is a solute in a 
liquid volume. 

24. An apparatus as claimed in claim 23, characterized 
in that said exchanger is a dialyzer comprising a semi-permeable 
membrane, where said solute and exchange liquid passes on one 
side of said membrane and said exchange fluid is a dialysis 
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liquid passing on the other side of said membrane resulting in a 
dialysate liquid exiting the dialyzer. 

25. An apparatus as claimed in claim 24, characterized 

by 

5 a flow rate measuring means for measuring the flow rate Q d of 
said dialysate liquid flow, and a third calculation means for 
calculating the accumulated mass of said solute in said 
dialysate liquid by integrating the product of said flow rate Q d 
and said solute concentration c d over time. 

10 26 • An apparatus as claimed in claim 24 or 25, 

characterized in that said dialysis liquid entering the dialyzer 
comprises a non-zero initial concentration of said composition, 
and that said measuring apparatus measures a concentration 
difference across the dialyzer. 

15 27. An apparatus as claimed in claim 24, 25 or 26, 

characterized in that said first calculation means is adapted 
for subtracting a compensation term (G/Q d ) from said solute 
concentration (c d ) to obtain a compensated concentration, taking 
the logarithm of said compensated concentration and fitting a 

20 straight line to said logarithm of said compensated 
concentration, where said compensation term compensates for 
generation (G) of said solute. 

28. An apparatus as claimed in any one of claims 22 to 

27, characterized in that said parameter of said approximation 
25 curve is a slope of said substantially straight line. 

29. An apparatus as claimed in any one of claims 22 to 

28, characterized in that said first calculation means is 
adapted to exclude data obtained during an initiation period, of 
for example 60 minutes. 

30. An apparatus as claimed in any one of claims 24 to 

29, characterized in that said first calculation means is 
adapted to adjust the time scale, when the flow of dialysis 
liquid is interrupted for a first time period, by replacing said 
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first time period with a replacement time period which is 
shorter than said first time period. 

31. An apparatus as claimed in any one of claims 24 to 
30, characterized in that said second calculation means 

5 calculates said parameter as the slope (K/V) with respect to 
time (t) of equation: 

ln(c d -G/Q d ) = ln(c 0 -G/Q d ) - Kt/V (I) 

where 

It c d = dialysate concentration at time t 

G = solute generation 
Qa = dialysate flow 

c 0 = dialysate concentration at time zero 
K/V = relative dialysis efficiency 
15 t = time from time zero. 

32. An apparatus as claimed in claim 31, characterized 
in that said second calculation means calculates a momentary 
mass (mi) according to the equation: 

mi = (Q d - Cdi) / (K/V)! (ii) 
20 where (K/V)= is determined according to equation (I), and in 
that the momentary relative efficiency (K/V} 2 at any time is 
determined according to equation: 

(K/V) 2 = (Qa ' c d2 ) / m 2 (III) 

where 

25 m 2 = m x - (U 2 - Ui) + G (t 2 - ti) (IV) 

and 

c d i = concentration at time ti 
c d2 = concentration at time t 2 
Ui = accumulated mass to time ti 
30 U 2 = accumulated mass to time t 2 . 

33. An apparatus as claimed in claim 32, characterized 
in that said second calculation means is adapted to integrate 
the momentary relative efficiency (K/V) over time to give an 
estimate of the total dialysis dose Kt/V. 
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34. An apparatus as claimed in claim 22, characterized 
in that said second calculation means is adapted to fit said 
approximation curve by calculating a line which passes through 
the greatest number of points in the logarithm of the 
concentration curve, possibly compensated. 

35. An apparatus as claimed in claim 34, characterized 
in that said second calculation means is adapted to use all 
values of concentration which are within a certain limit from 
said line for calculating the momentary masses m„, which are 
then used for calculating the Initial masses m 0n , and in that 
the calculated initial masses m 0n are used for estimating the 
actual initial mass m 0 , for instance by taking the median or 
mean value of the initial masses m 0n . . 

36. An apparatus as claimed in claim 32 or 33, 
characterized by means for determining a momentary mass m : in 
another way, for example by analysing a blood sample or by 
equilibration of dialysate with blood and determining the actual 
concentration of the composition and by estimating or measuring 
the distribution volume of the composition. 

37. An apparatus as claimed in anyone of claims 24 - 34, 
characterized by means for estimating the distribution volume 
(V) of the composition fluid, for example by Watson's formula, 
and that the concentration cf the composition in the fluid is 
determined by dividing the calculated mass with the volume. 

38. An apparatus as claimed in anyone of claims 24 - 36, 
characterized by means for measuring the concentration of the 
composition in the fluid and that the distribution volume is 
determined by dividing the calculated mass and the 
concentration. 

39. An apparatus as claimed in anyone of claims 24 - 36, 
characterized by means for measuring the concentration of the 
composition in the fluid comprising means for introducing a 
disturbance into the dialyser and means for measuring the 
resulting effect in the effluent dialysate and means fcr 
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calculating the effective clearance of the dialyser from this 
resulting measurements, and means for calculating the plasma 
water concentration of said substance by the formula: 
Cpw = Q d x c d / K e 

where 

c pw =plasma water concentration of urea at initiation of dialysis 
Q d = effluent dialysate flow rate 

c d = concentration of urea as extrapolated to the initiation 
K e = effective clearance of the dialyser for urea 
and means for determining the distribution volume v of said 
substance by the formula: 
V = m 0 / c pw 

where m 0 is determined according to any of claims 24-38. 

40. An apparatus as claimed in anyone of claims 22 - 33, 
characterized by means for determining a deviation of the 
concentration curve from the approximation curve, after 
establishment of said approximation curve, and emitting an alarm 
at deviation over or under a predetermined threshold level. 

41. An apparatus as claimed in any one of claims 24 to 
30, characterized in that said second calculation means 
calculates said parameter as the slope [(K-UF)/V 0 ] of ln(c d -c k ) as 
a function of (V 0 /UF) In (V/V 0 ) in the equation: 



where 



ln(Cd-Ck) = ln(c 0 -c k )+ [ (K-UF) /V 0 j- [ (V 0 /UF) In (V/V 0 ) ] (V) 

c d = dialysate concentration at time t 

c k = G/ [Q d (1 -UF/K) ] 

G = solute generation 

Qd = dialysate flow 

c 0 = dialysate concentration at time zero 

K/V = relative dialysis efficiency 

UF = ultrafiltration per time 
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42. An apparatus as claimed in claim 41, characterized 
in that said second calculation means calculates a momentary 
mass (m x ) according to the equation: 

mi = (Q d " cai) / (K/Vh (VI) 

where 

(K/Vh = (K/V 0 ) / (1 - t x ■ OF/ V 0 ) (VII) 

with (K/V 0 ) determined according to equation (V) and UF/V 0 
estimated, and in that the momentary relative efficiency (K/V; 2 
at any time is determined according to equation: 

(K/V) 2 = (Q d - Cd2 ) / m 2 (VIII) 



where 
and 



m 2 = mi - (U 2 - Ui) + G (t 2 - t a ) (ix) 

Cdi = concentration at time ti 
Cd2 = concentration at time t 2 
U 2 = accumulated mass to time ti 
U 2 = accumulated mass to time t 2 . 
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